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1. Executive Summary 

This report summarizes the research carried out as part of the project "An Optimization 
Model for the U. S. Air-Traffic System (NASA# NAG- 1-520). The primary research objec- 
tive was to establish a systematic approach for monitoring U. S. air traffic in the context of 
system-wide planning and control. Towards this end, a network optimization model with non- 
linear objectives was chosen as the central element in the planning/control system. The net- 
work representation was selected because: (1) it provides a comprehensive structure for dep- 
icting essential aspects of the air traffic system, (2) it can be solved efficiently for large scale 
problems, and (3) the design can be easily communicated to non-technical users through com- 
puter graphics. 

Briefly, the network planning models consider the flow of traffic through a graph as the 
basic structure. Nodes depict locations and time periods for either individual planes or for 
aggregated groups of airplanes. Arcs define variables as actual airplanes flying through space 
or as delays across time periods. As such, a special case of the network can be used to model 
the so called flow control problem 2 . 

Due the large number of interacting variables and the difficulty in subdividing the prob- 
lem into relatively independent subproblems, we designed an integrated model which will 
depict the entire high level ( above 29000 feet) jet route system for the 48 contiguous states in 
the U. S. While the resulting model is large by today’s standards, we felt that computer tech- 
nology (hardware and software) is improving rapidly and a practical full-scale system should 
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• become feasible over the next few years. 

As a first step in demonstrating the concept’s feasibility, we gathered data from the Indi- 
anapolis control sector and built a nonlinear risk/cost model for this airspace. The nonlinear 
network program — NLPNETG- was employed in solving the resulting test cases 1 . This 
optimization program uses the Truncated-Newton method (quadratic approximation) 3 for 
determining the search direction at each iteration in the nonlinear algorithm. 

It was shown that aircraft could be re-routed in an "optimal" fashion whenever traffic 
congestion increased beyond an acceptable level, as measured by the nonlinear risk function. 
An efficient cost/risk frontier can be traced out by altering the relative weights between cost 
and risk. This curve can be employed within the context of strategic planning, e.g. adding air- 
port resources, or as an operational decision tool for assisting air-traffic controllers. 

Computational costs for this exercise were quite reasonable and the nonlinear optimiza- 
tion problems could be solved on an minicomputer - VAX 11/750 —at Princeton University. 
In addition, we transported the optimization system to a CRAY supercomputer at Boeing 
Computer Services, Seattle, Washington. Taking advantage of the CRAY vector architecture 
required a substantial effort; see reference 6 for a description of the steps that were taken. This 
analysis demonstrated that the Trucated-Newton algorithm can be tuned for a vector com- 
puter and that substantial air-traffic control problems could be solved in real-time. 

The technical details of the project are contained in three research reports. The report 
titles are listed below: 

(1) "Nonlinear Network Programming on Vector Supercomputers," Report EES-85-13, 
Princeton University, submitted to Operations Research , 1986. 

(2) "Real-Time Operational Planning for the U.S. Air-Traffic System", Report EES-86-5, 
Princeton University, submitted to Applied Numerical Mathematics , 1986. 

(3) "Integrated Risk/Cost Planning Models for the U. S. Air Traffic System," Report EES- 
85-9, Princeton University, submitted to Management Science , 1985. 
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These reports are self-contained, describe the research work performed during the project, 
and are provided in the Appendix. 

2. Future Research 

We believe that future research involving integrated air-traffic modeling should be 
directed along several avenues. First, the size of the network test problems should be 
expanded by including multiple control sectors and a longer time horizon -- more periods. As 
the problems grow in size, every attempt must be made to keep the execution time within a 
reasonable level. The runtime issues are especially pertinent in the case of the operational 
planning model. Remember that the final goal is to develop a practical decision support sys- 
tem for air traffic planning. 

Second, an interface should be designed so that non-technical personnel will be able to 
effectively use the planning systems. It will be essential to understand the type and the form of 
information which is needed by the air traffic controller, among others. An interactive pro- 
cedure which employs graphics seems to be the most likely format. Perhaps, color and inten- 
sity could be used to identify the critical elements in the control domain. Also, the optimiza- 
tion procedure needs to be flexible enough to be able to handle alternative objective func- 
tions, priorities and possible intervention by the users. Criteria other than utilitarian can be 
achieved by adjusting the objective function in the network models. 

Third, the risk objective function which measures issues such as congestion, proximity, 
and workload, must be refined and tested with empirical findings. There is much to be done 
in this arena. For example, the analysis must be coordinated with the micro-level studies car- 
ried out by Odoni 6 and others. The field of risk analysis is only beginning to develop general 
principles. As an approximation, one could use a congestion function which depends upon 
geometry, time of day, weather, and other factors to evaluate the degree of acceptability for an 
airspace as compared to delaying or re-routing planes anticipated to cross the region of 
interest. See the second and third papers listed in the Appendix for further details. 
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In summary, the National Airspace Plan 4 promises to greatly enhance the quality and 
the quantity of information available on the state of the U.S. air traffic system. Never die less, 
much effort will be needed to decide what information is critical to a safe and efficient air- 
traffic environment. In this regard, planning models will assist in important task of identifying 
critical informational flows. Both the strategic and the operational network models proposed 
in this project are designed for this task. 
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This appendix contains the computer program used for generating the 
Air-Traffic test cases in the project. 

c*** ********************************************************************* 
£*********************^************************************************** 

‘ £ -k * * * * 

C** PROGRAM FAANET *** 

Q * * * * * 

Q'k-k-k-k-k’k’k-k-k-k-k-k'k'kic'k'k'k'k'k'k'k'k'k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k'k-k-k-k’k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k'k 

c ************************************************************************ 

C 

c PURPOSE : 

c This program reads the data for the airports, flights and 

c fuel burn model and generates the network model for input to the 

c nonlinear network optimizer NLPNETG 

c It compiles with F77 on the VAX 11/750 under UNIX 5.2 

c 


c SUBROUTINES : 

c 

c 

c 

c 

c 

c 

c 

c INPUT FILES : 

c 

c 

c 

c 

C OUTPUT FILES: 

c 
c 
c 
c 

c* 
c 

c- 
c 

c . 
c 
c 
c 
c 
c . 
c 


INPUT : Reads in the data for Airports/Flights/Fuel burn 
GENER : Generates the network topology and identifies risk 

BODY : Writes on the output file the generated network 

TAIL : Writes on the output file the risk function 

LENGTH : Calculates the distance traveled by a flight 


Logical unit 7 
Logical unit 8 
Logical unit 9 


Logical unit 10 
Logical unit 15 
Logical unit 16 


Fuel burn model data 
Airports location data 
Flights strips data 


Nonlinear network model 

Warning/Error messages from the program 

Scratch file for storage of interractions 


a********************************************************************** 
-- DEFINITION OF VARIABLES 
. . CHARACTER Arrays 


airprt (nairpt) 
fights (nf light) 

REAL* 4 Arrays 


Four character airports code 
Six character flight id code 


c 

burn 

(nflight *2, 9) : 

Fuel burning rates at 9 legitimate 

cruise altitudes 

c 

cost 

(narc) : 

Cost of generated network arcs 


c 
c . 

. . INTEGER 

Arays 



c 

c 

iorig 

(nflight) : 

Pointer to origin airport for every flight 

c 

idest 

(nflight) : 

Pointer to destination airport 

it 

c 

if lght (nflight) : 

Pointer to fuel burn record 

ti 

c 

alti 

(nflight,2) : 

Entering/Exiting altitude 

ii 

c 

itime 

(nflight, 2) : 

Entering time ( hour/min) 

»i 

c 

ot ime 

(nflight, 2) : 

Exiting time ( hour/min) 

it 

c 

c 

hemi 

(nflight) : 

Hemi code 

ti 

c 

rate 

(nflight) : 

Burning rate indicator 

ti 

c 



- 1 Lbs/nmi 


c 



- 2 Lbs /hr 


c 

c 

lati 

(nairpt, 3) : 

Latitude of every airport 


c 

long 

(nairpt, 3) : 

Longtitude of every airport 


c 

elev 

(nairpt) : 

Elevation of every airport 


c 

c 

ii 

(narc) : 

From node of generated network arcs 


c 

j j 

(narc) : 

To node of generated network arcs 
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c 

c istart(9) : Pointer to array ITAG for arcs contributing 

c to the risk function at different altitudes 

c itag (narc ) : Thread of arcs contributing to risk 

c 

c. . . INTEGER Variables 
c 

c nf light : Number of flights 

c itper : Number of time periods considered for flights delay 

c time : Length of each delay time period 

c inter : Number of time intervals considered for risk analysis 

c 

c * ********************************************************************** 
£*************************************************★********************* 

c 

c 

program faanet 
c 

c 

c 

include 7 - . /data/comdat 7 
c 

open (7) 
open (8) 
open (9) 
open (10) 
open (15) 
open (16) 
print *, 7 [H [J 7 
c 

c Read Input data 

c 

write (15,900) 
call input 
c 

c Generate and write out the Network topology 

c 

write (15, 910) 
call gener 
call body 
c 

c Determine and write risk interractions over target sector 

c 

write (15,920) 
call risk 
call tail 

c 

print *, 7 [H [J #3 #4 D 0 N E ! [2B 7 

c 

c FORMAT Statements 

c 

900 format </,lx, 7 *** READ INPUT DATA ... 7 ,/) 

910 format (/,lx, 7 *** GENERATE NETWORK TOPOLOGY ... 7 ,/) 

920 format (/,lx, 7 *** DETERMINE FLIGHT INTERRACTIONS 

c 

c 

close (7) 
close (8) 
close (9) 
close (15) 
close (16) 
c 

stop 

end 

c 

c *** ******************************************************************** 
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non oooo ooo oooonoo 


c 

subroutine input 

*********************************************************************** 
PURPOSE : 

This subroutine reads the data from the three 
input files as described in the main program . 

include ' . . /data/comdat ' 
character*6 aflt 
character*4 aorig, adest 

Prompt the user for control parameters 

print [ 1 ; 4 ; 7m #3 #4 #6 FAA NETWORK GENERATOR [Om [3B' 

print *, 'Number of time periods considered for flight delay ?' 
read (5,*) itper 

print Length of flight delay time interval ?' 
read (5,*) time 

print *, 'Number of time periods considered for risk analysis ?' 
read (5,*) inter 

print *, 'Value of printing flag ?' 

, read (5,*) mprint 

print *, 'Header card for generated network ( at most 72 char) ?' 
read (5,890) title 

print [ 3B [4C [l;4;5;7m #3 #4 #6 R U N N I N G [0m' 

Read airports code and location 

i = 1 

10 read (8, 900, end=20) airprt(i), lati (i, 1) , lati (i, 2) , lati(i,3), 

* long(i,l), long(i,2) ,long(i,3) ,elev (i) 
i = i + 1 

go to 10 

Read fuel burn data 

. 20 nairpt = i - 1 
i = 1 

c 

30 read (7, 910,end=40) fights (i) ,hemi (i) , rate (i) , (burn (i, j) , j=l, 9) 
i = i + 1 
go to 30 
c 

c . Read flights data and set pointers to airport and fuel-burn recors 

c 

40 nfbrec = i - 1 
i = 1 
c 

50 read ( 9, 920, end=100) aflt, aorig, adest , ialtl, ialt2 , 

* (itime.(i, j) , j=l,2) , (otime(i,j), j=l,2) 
c 

alti (i, 1) = ialtl * 100 
alti (i, 2) = ialt2 * 100 
c 

iflagl = 1 
if lag2 = 1 
c 
c 

do 60 j = 1, nairpt 

if (aorig. eq.airprt ( j) ) then 
iorig(i) = j 
iflagl = 0 

end if 
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if (adest .eq.airprt ( j) ) then 
idest(i) = j 
iflag2 = 0 

end if 

if (iflagl.eq.O .and.iflag2 .eq.O) go to 70 

60 continue 
c 

if (if lagl . ne . 0) write (15,930) aorig , aflt 
if (if lag2 . ne . 0) write (15,935) adest , aflt 

c 

70 do 80 j = 1, nfbrec 

if (af It -eq. fights (j) ) then 
iflght(i) = j 
go to 90 

end if 

80 continue 

c 

write (15,940) aflt 

c 

90 i = i + 1 

go to 50 

100 nflight = i - 1 
c 

if (mprint . le . 3) go to 999 
write (15,950) 
write (15,960) 

write (15,970) (i, iorig (i) , idest (i) , 

* (alti (i, j) , itime (i, j) ,otime (i, j) , j=l,2) , i=l, nflight) 
write (15, 980) 

write (15, 990) (airprt (i) , lati (i, 1) , lati (i, 2) , lati (i, 3) , 

* long (i, 1) , long (i, 2) , long (i, 3) , elev (i) , i=l, nairpt) 
write (15,992) 

write (15,994) (i, fights (i) , if lght (i) , hemi (i) , rate (i) , 

* burn (i, 1) , i=l, nfbrec) 

c 

c FORMAT Statements '■ 


890 


format (a72) 

900 


format 

(lx 

910 


format 

(lx 

920 


format 

(lx, 

930 


format 

(lx 


* 

r 

of 

935 


format 

(lx 


•k 


a4 

940 


format 

(lx 


i k 


lx 

950 


format 

(lx 

960 


format 

</, 


k 


,2 

970 


format 

(lx 

980 


format 

</, 

990 


format 

(2x 

992 


format 

(lx 


k 


5x 

994 


format 

(lx 

999 


return 




end 



c 

c 


block data 


c 
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c 

include ' . . /data/comdat ' 
c 

data ihO / 29000 , 33000 , 31000 , 41000, 45000/ 
data ihl 731000,35000,39000,43000/ 
data nfaa/0/ 

‘ c 

end 
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c 

subroutine gener 
c 

c 


c PURPOSE 

c 

c 

c 

c 

c 

c 

c 

c INPUT 

c 

c 

c 

c 

c OUTPUT 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


This subroutine generates the network topology for 
the Flight/Airport information given for one control sector. 

It generates alternative cruise altitudes for every flight 
and determines associated flight costs as well as delay costs . 

For every generated route (time/altitude) its contribution 
to the risk function is determined . 


The arrays read from the database files by routine INPUT 
ihO(nflight) , legitimate cruise altitudes for hemi-code = 0 
ihl(nflight) , legitimate cruise altitudes for hemi-code = 1 


ii (narc) 
jj (narc) 
cost (narc) 
ipntf 1 (narc) 
iper (narc) 


from node of generated arc 
to node of generated arc 
cost of generated arc 
flight number of this arc 


, time period of this arc 


istart (inter*9) , pointer to ITAG array 


itag (narc) 
narc 


thread of lfights contributing to risk 
total number of generated arcs 


include ' . . /data/comdat ' 


narc =* 0 
i =0 
j - 0 


c== Consider origin nodes for all time periods. ===== 

c 
c 
c 


do 1000 i = 1, itper 

do 500 j = 1, nflight 


iinod = (i-1) *nf light + j 

jjl = i*nflight + j 

jj2 = (itper+i-1) *nf light + j 

itin = itime ( j, 1) *60 + itime(j,2) 
itout = otime ( j, 1) *60 + otime(j,2) 


c Takeoff delay arcs 

c 


(i-1) *time 
(i-1) *time 


ipnt = iflght (j) +1 

if (rate (ipnt) .ne. 2) write (15/900) fights (ipnt) 
costd = time * burn (ipnt/ 1) / 60.0 


c En route arcs 

c 


ipnt = iflght (j) 
ihcode= hemi (ipnt) 
ialtO = 0 

ialtl = alti (j/1) 
ialt2 = 0 
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c 

c 


determine alternative cruise altitudes for this flight 


c 

c 

c 


c 


20 

30 

c 


40 

c 

45 


c 


c 



70 

80 

c 


90 

c 

95 


c 


c 


if (ihcode .eq. 1) go to 50 

alternatve routes for Hemi code 0 flights 

if (alti(j,l) .eq. alti(j,2)) go to 30 
ialtl = ( alti(j,l) + alti(j,2)) / 2 

do 20 k = 1, 4 

if (ialtl . It . ihO (k+1) . and. ialtl . gt . ihO (k) ) then 
ialtl = ihO (k+1) 
go to 30 

end if 
continue 
continue 

do 40 k = l f 5 

if (ialtl .eq. ihO (k) ) go to 45 
continue 

fburnl = burn (ipnt,2*k-l) 

if (fburnl .eq. 0 . 0) write (15, 910) fights (ipnt) 

if (k.gt.l) then 

ialtO = ihO(k-l) 

fburn0= burn (ipnt,2*k-3) 

if (fburnO .eq. 0 . 0) fburnO = fburnl 

end if 

if (k.lt.5) then 

ialt2 = ihO (k+1) 

fburn2= burn (ipnt,2*k+l) 

if (fburn2 .eq. 0 . 0) fburn2 = fburnl 

end if 
go to 100 


alternatve routes for Hemi code 1 flights 
continue 

if (alti ( j , 1) .eq. alti ( j , 2) ) go to 80 
ialtl = ( alti ( j, 1) + alti(j,2)) / 2 

do 70 k = 1, 3 

if (ialtl . It . ihl (k+1) .and. ialtl .gt . ihl (k) ) then 
ialtl = ihl (k+1) 
go to 80 

endif 

continue 

continue 

do 90 k = 1, 4 

if (ialtl .eq. ihl (k) ) go to 95 

continue 

fburnl = burn (ipnt , 2*k) 

if (fburnl .eq. 0 . 0) write (15, 910) fights (ipnt ) 

if (k.gt.l) then 

ialtO = ihl(k-l) 

fburn0= burn (ipnt,2*k-2) 

if (fburnO .eq. 0 . 0) fburnO = fburnl 

end if 
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o o 


c 

c . . . 
c 

100 

c 


c 

c— 

c 

c : 

c 


c 

c 

c 

c 


c 

c 


c 

500 

1000 

c 

c== 


if (Jc.lt. 4) then 

,ialt2 = ihl(k+l) 
fburn2= burn (ipnt,2*k+2) 
if (fburn2 .eq. 0 . 0) fburn2 = fburnl 

end if 

determine the cost of alternative altitudes 
continue 

call length (j ,dist) 

costO = fburnO * dist 

costl = fburnl * dist 

cost2 = fburn2 * dist 


Store generated arcs. (For the last time period no delays possible) 


> Save Delay arcs 

if (i.ne.itper) then 

narc = narc + 1 
ii (narc) = iinod 
jj (narc) = jjl 
cost (narc) = costd 
ipntfl(narc) = j 
iper (narc) = i 

end if 

Save flight arcs on different altitudes 

PRIMARY-1 

if (ialtO.gt.O) then 

narc = narc + 1 
ii (narc) = iinod 
jj (narc) = jj2 
cost (narc) = costO 
ipntfl(narc) = j 
iper (narc) = i 

end if 


PRIMARY 

narc = narc + 1 
ii (narc) = iinod 
jj (narc) = jj2 
cost (narc) = costl 
ipntfl(narc) = j 
iper (narc) = i 


PRIMARY+1 

if (ialt2.gt.O) then 

narc = narc + 1 
ii (narc) = iinod 
jj (narc) = jj2 
cost (narc) = cost2 
ipntfl(narc) = j 
iper (narc) = i 

end if 

continue 

continue 


Consider destination nodes for all time periods 
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o o 


c 

itor = itper * nflight 
itort2 = itor * 2.0 

c 

do 2000 i = 1, itper 

. c 

do 1500 j = 1, nflight 

c 

iinod = itor + (i-1) *nf light + j 
jjl = itor + i *nflight + j 
j j2 = itort2 + j 

c 

c Landing delay arcs . . . 

c 

ipnt = iflght (j) +1 

if ( rate (ipnt) .ne . 2) write (15,900) flghts(ipnt) 
costd= time * burn (ipnt, 1) / 60.0 

c 

c Store generated arcs. (For the last time period no delays possible) 

c 

if (i.ne. itper) then 

narc = narc + 1 
ii(narc) = iinod 
j j (narc) = jjl 
cost (narc) = costd 
ipntfl(narc) = j 
iper (narc) = i 

end if 

c 

narc = narc + 1 
ii(narc) = iinod 
j j (narc) = j j2 
cost (narc) = 0.0 
ipntfl(narc) = j 
iper (narc) = i 
c 

1500 continue 

2000 continue 
c 

c FORMAT Statements 

c 

900 format (lx, 'Check fuel burn data for flight ',a6) 

910 format (lx, 'Fuel burn data for flight ',a6,' is not available') 


return 

end 


c 

c 


subroutine length (iflt, dist) 


c 

c **** ******************************************************************* 
c 

c PURPOSE : 

c Returns the distance between two airports - origin and 

c destination - of a flight. 

c NOTE : The distance is measured ignoring difference in 

c the elevation of the origin-destination airports and the cruise 

c altitude of the flight. The earth curvature is taken into account. 

c 


c INPUT : 

c iflt 

c airprt (nairpt) 

c lati (nairpt, 3) 

c long (nairpt, 3) 


Flight number indicator 
List of airport codes 
Latitude of all airports 
Logitude of all airports 
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c iorig (nflght) : Pointer to origin airport for flight IFLT 

c iciest (nflght) : Pointer to destination airport for flight IFLT 

c 

c degrad : PI/180 - converts degrees into radians 

c radmil : 24844mls /2*PI*1.151 -converts rad into naut . miles 

c 

include ' . . /data/comdat ' 
real latl, lat2 
c 

c Parameter initialization 

c 

parameter ( degrad = 0.017453 , radmil = 3439.4 ) 
ipo = iorig (iflt) 
ipd = idest (iflt) 
c 

c 

c Calculate latitudes of Origin/Destination airports in radians 

c Compute the linear distance between the airports as : 

c (GREAT CIRCLE DISTANCE) * (EARTH CIRCUITY FACTOR) 

c 

c 

xl = lati(ipo,l) + lati(ipo,2) / 60.0 + lati(ipo,3) / 3600.0 

yl = long(ipo,l) + long(ipo,2) / 60.0 + long(ipo,3) / 3600.0 

x2 = lati(ipd, 1) + lati(ipd, 2) / 60.0 + lati(ipd,3) / 3600.0 

y2 = long(ipd^l) + long(ipd,2) / 60.0 + long(ipd,3) / 3600.0 

c 

latl = degrad * yl 
lat2 = degrad * y2 
c 

dist = acos (sin (latl) *sin (lat2) + cos (latl) *cos (lat2 ) * 

1 cos (degrad* (xl-x2) ) ) *radmil 

c 

if (mprint .ge . 3) write (15,900) iflt, dist 

c 

c FORMAT Statements 

c 

900 format (lx, 'Length of flight ',i3,' is f ,e25.12) 

c 

c 

return 

end 
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C 

subroutine body 
c 

c * * * * ******************************************************************** 

, c 

c PURPOSE : 

c This subroutine writes on the output file the body 

c of the generated network, 

c 

c INPUT : 

c The arrays generated by routine GENER 

c 

c 

c 

include ' . . /data/comdat ' 

c 

c Initialization 

c 

character*l a, s , d 
character*2 bd 

parameter (i0=0, il=l, i2=2, r0=0, rl=l) 
data a/lha/, s/lhs/, d/lhd/, bd/2hbd/ 

. c 

c Check the size of generated network 

c 

nodes = 2*nf light*itper + nflight 

na = 3*nf light*itper + 2*nf light* (itper-1) + nflight *itper 
nn = jj (narc) 
c 

if (na.ne.narc) write (15, 900) na, narc 
if (nn .ne .nodes) write (15,910) nodes, nn 
c 

c Write header card 

c 

write (10,920) title 
write (10,930) nn, narc, i2, il 
c 

c Write arc data 

c 

write (10,940) a 
do 100 i = 1, narc 

write (10,950) bd, iO, ii(i), j j (i) , rO , rl, rO, cost ( i) , rO , rO 

100 continue 
c 

c Write supply data 

c 

write (10,940) s 
do 200 i = 1, nflight 

write (10,970) i , rl 

200 continue 

c 

c Write demand data 

c 

write (10,940) d 
itemp = 2*nf light*itper 
do 300 i = 1, nflight 

nod = itemp + i 
write (10,970) nod, rl 

300 continue 

c 

c FORMAT Statements 

c 

900 format ( lx, ' WARNING : Model should have ' , i5,' arcs. Only ',15, 

* ' were generated. ', /, 

* llx, ' Data for some flights may be missing. Check the databases',/) 
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910 


format (lx, 'ERROR : Model should have nodes. Only ',i3, 

* ' were generated.',/, 

*9x, ' Possible program error.',/) 


920 

format 

(lx,a72) 

930. 

format 

(4i5) 

940 

format 

(al, t72, ' ') 

950 

format- 

(3x,a2, 3x,i2,2i5,3f 10,1, 3f 10.1) 

970 

format 

(10x,i5,5x,fl0.0,t72, ' ') . 


return 

end 



c 

c************************************************************************* 

c 

subroutine tail 
c 

c************************************************************************* 

c 

c PURPOSE : 

c This subroutine writes on the output file the components 

c of the risk function (ie. nonlinear and non-separable function) . 

c For the special case when only a single flight is crossing 

c the sector the risk function is separable. 

c 

c INPUT : The arrays generated by routine RISK 

c 

c 

include ' . . /data/comdat ' 

c 

character*3 f , e 
parameter (zero=0.0, two-2.0) 
data f/3hfaa/, e/3hend/ 
imax = - 1 
imin = 99999 
rmxcst= 0.0 
c 

c > Determine the max cost to use as scaling factor. 

c 

do 5 i = 1, narc 

costi ~ cost (i) 

if (rmxcst .It .costi ) rmxcst - costi 

5 continue 

rmxcst = rmxcst / 10.0 

c 

c=== Determine the planning interval 

c 

do 10 i = 1, nf light 

itmax = otime (i, 1) *60 + otime(i,2r 
itmin = itime (i, 1) *60 + itime (i, 2) 
if (itmax. gt .imax) imax = itmax 
if (itmin.lt .imin) imin = itmin 

10 continue 

c 

c > Length of interval used for risk evaluation 

c 

rintl = real (imax-imin) / real (inter) 
write (15,900) imin, imax, rintl 

c 

c=== Consider all flights at the same cruise altitude 

c 

write (10,910) f 
do 1000 i = 1, 9 
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c 

ipi = istart (i) 

if (ipi .eq. 0 .or . itag (ipi) .eq. 0) go to 1000 
c 

100, if (ipi.eq.O) go to 1000 

ip j = ipi 

200 ip j = itag (ipj) 

if (ipj.eq.0) go to 500 

c 

iarci = ipntfl (ipi) 
iarcj = ipntfl (ipj) 

c 

if {iarci. eq. iarcj) go to 200 

c 

iaddi = time * (iper(ipi) - 1) 
iaddj = time * (iper(ipj) - 1) 
isti = itime (iarci, 1) *60 + itime (iarci, 2) 
iendi = otime (iarci, 1) *60 + otime (iarci, 2) 
istj = itime (iarcj, 1) *60 + itime (iarcj, 2) 
iendj = otime (iarcj , 1) *60 + otime (iarcj, 2) 
c 

ist = maxO (isti , istj ) 
iend = minO (iendi, iendj) 
if (ist .ge. iend) go to 200 


c 

c > ALPHA is proportional to the overlapping time in sector 

c 

rtemp = iend - ist 
alpha = rmxcst * (rtemp / rintl) 
write (16,920) ipi, ipj, alpha, two, zero 
go to 200 
c 

500 ipi = itag (ipi) 

go to 100 

c 

1000 continue 


+ .iaddi 
+ iaddi 
+ iaddj 
+ iaddj 


c 

write (16,920) narc, narc, zero, zero, zero 
close (10) 

call system ('sort +10 -15 < fort. 16 » fort. 10') 
c write (10,910) e 

if (mprint .ge.2) write (15, 930) 

* (istart (i) , i=l, 9) , (itag (i) , i=l, narc) 

c 

c FORMATS 

900 format ( lx, ' Planning interval: From ',i3,'min to ',i4, 

* lx, 'min. Length of risk interval f 6 . 2) 

910 format (a3, t25, ' 1 . 0 ' , t72 , ' ') 

920 format (10x, 2i5, 3f 10 . 4) 

930 format (lx, ' ISTART (i) ' , 9 (lx,i5) , /, lx, ' ITAG (i) ' , / , 100 ( 10 ( lx, i5) , / ) ) 

c 

return 

end 
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o o 


*********************************************************************** 

subroutine risk (num, ist) 
c 

£*********************************************************************** 

C 

c PURPOSE : 

c This subrutine generates the arrays for the nonlinear 

c part of the network model. For a given flight/time it identifies the 

c altitude and time interval when a conflict may arise and creates the 

c thread that links all interracting flights 

c 

c INPUTS : 

c 
c 
c 

c OUTPUTS 

c 
c 
c 
c 

C' 

include 9 . . /data/comdat ' 
c 

c=== Detemine the flights that cross the target sector at the same altitude ===== 
c 

c > Initialize storage arrays 

c 

do 12 i = 1, 9 

istart (i) = 0 

12 continue 

do 15 i = 1, narc 

itag (i) = 0 

15 continue 


num , number of risk periods while flight is in target sector 
ist , pointer to the flight altitude while in target sector 

istart, pointer to the starting location of the thread of 
interracting flights 
itag , thread of interracting flights 


do 500 i = 1, itper 
do 200 j = l f nflight 

c 

ipnt = iflght (j) 
ihcode= hemi (ipnt) 
ialtO = 0 

ialtl = alti (j,l) 
ialt2 = 0 
c 

iinod = (i-1) *nflight + j 
j j 2 = (itper+i-1) *nf light + j 

c 

c > Find all arcs for this flight 

c 

icount = 0 

18 icount = icount + 1 

if (ii (icount) .eq. iinod. and. j j (icount) .eq. j j 2 ) goto 19 
if (icount . It .narc) go to 18 
write (15,910) fights (j), i 
go to 200 
c 

19 continue 
c 

c > Determine the flight altitude 

c 

if (ihcode .eq. 1) go to 50 
c 

c... Altitude for Hemi code 0 flights 
c 


- 14 - 


c- 


if (alti ( j, 1) .eq.alti ( j, 2 ) ) go to 30 
ialtl = (alti(j,l) + alti(j,2)) / 2 

do 20 k = 1, 4 

if (ialtl . It . ihO (k+1) . and. ialtl .gt . ihO (k) ) then 
ialtl - ihO (k+1) 
go to 30 

end if 

20 continue 

30 continue 

do 40 k = 1, 5 

if (ialtl .eq. ihO (k) ) go to 45 
40 continue 

45 continue 

--> Save all the arcs corresponding to this flight at all altitudes 

— PRIMARY - 1 

if (k .ne . 1 .and. ii (icount) . eq. iinod. and. j j (icount ) . eq. j j2) then 
kml = k - 1 

call save (kml, icount) 
icount = icount + 1 

end if 

— PRIMARY 

if (ii (icount) . eq. iinod. and . j j (icount ) .eq. j j 2 ) then 
call save (k, icount) 
icount = icount + 1 

end if 

-- PRIMARY + 1 

if (k.ne.5.and.ii (icount) . eq. iinod. and . j j ( icount ) . eq . j j2 ) then 
kpl = k + 1 

call save (kpl, icount) 

end if 

go to 200 

. . Altitude for Hemi code 1 flights 

50 continue 

if (alti ( j, 1) .eq.alti ( j , 2) ) go to 80 
ialtl = (alti ( j , 1) + alti(j,2)) / 2 

do 70 k = 1, 3 

if (ialtl . It . ihl (k+1) . and. ialtl . gt . ihl (k) ) then 
ialtl = ihl (k+1) 
go to 80 

end if 

70 continue 

80 continue 

do 90 k = 1, 4 

if (ialtl .eq. ihl (k) ) go to 95 
90 continue 

95 continue 

— > Save all the arcs corresponding to this flight at all altitudes 

— PRIMARY - 1 

if (k .ne . 1 .and. ii (icount ) .eq. iinod. and. j j (icount ) .eq. j j 2 ) then 
kml = (k - 1) + 5 
call save (kml, icount) 
icount = icount + 1 

end if 

- PRIMARY 

if (ii(icount) . eq. iinod . and. j j (icount ) .eq. j j 2 ) then 
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kp5 = k + 5 

call save (kp5, icount) 

icount = icount + 1 

end if 

C— t PRIMARY + 1 

if (k.ne . 4 .and. ii (icount) . eq. iinod. and . j j (icount ) .eq. jj2) then 
kpl = (k + 1) + 5 
call save (kpl, icount) 

end if 
c 

200 continue 

500 continue 

c 

c FORMATS 

905 format (lx, 'Flight ',a6,' at time period', i3, 

* ' is in target sector. Time in:',i3,' Time out:',i3) 

910 format (lx,'*** WARNING *** No arcs were generated for flight 

* a6,' at time period ',i5) 
c 

return 

end 


c 

subroutine save (ipos, iarcn) 
c 

c * ********************************************************************** 

c 

c PURPOSE : 

c This subroutine stores the specified flight 

c in the array of all flights that are crossing the target sector 

c simultaneoulsy and at the same level 

c 

c INPUTS : 

c iarcn , arc number for flight (time/alittude) considered 

c ipos , pointer to array ISTART indicating altitude 

c— OUTPUTS : 

c istart , pointer to thread storing all flights at this level 

c itag , thread of all flights crossing target sector at the 

c same interval 

c 

include ' . . /data/comdat ' 
c 

if (mprint .ge . 3) write (15,900) iarcn, ipos 

c 

c Is this the first flight at this altitude ? 

c 

ip - istart (ipos) 
if (ip.eq.0) then 

istart (ipos) = iarcn 
go to 9999 

end if 

ipml = ip 

c - 

c If more flights at this altitude exist find an empty spot to store this 

c 

100 ip = itag (ipml) 
if (ip.eq.0) then 

itag (ipml) = iarcn 
go to 9999 

end if 

ipml = ip 
go to 100 
c 

9999 if (mprint .ge . 3) then 

write (15,910) (istart (i) , i=l, 9) 
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write (15,920) (itag(i), i=l,narc) 

end if 
c 

c FORMATS 

906 format (lx, 'Save arc No. ',i5,' in altitude position ',i3) 
910 format (lx, ' ISTART (i) :',9(lx,i5)) 

920 format (lx,'ITAG(i) : ' , 100 (10 (lx, i5) , /) ) 

c 

return 

end 
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10 


character*28 char(l) 
character*48 char2(l), char3(l) 
read (7 , 100 , end=500 ) char(l) 
read (8,110) char2(l) 
read (9,110) char3(l) 

write (10,200) char (1) (2 : 8) , char (1) (9 : 12) , char ( 1 ) ( 13 : 14 ) 

, char (1) (15:15) , char (1) (16 : 16) , char (1) (17:29) 

* , char2 ( 1 ) ( 1 : 13) , char2 (1) (13 :25) , char2 (1) (25 : 37 ) , char2 ( 1 ) (37:49) 

* , char3 (1) ( 1 : 13 ) , char3 (1 ) (13:25) ,char3 (1) (25 : 37 ) , char3 ( 1) (37:49) 

100 format (a28 ) - 

110 format (a48) 

20 0 format (lx, a6, 3x, 2x, a3, 3x, a2, 2 (lx, al) , 9 (al2) ) 
go to 10 
c 

500 print *, 'Finished' 

stop 
end 
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10 


► * 
* 

100 

200 

c 

500 


character*21 char(l) 
read (9, 100,end=500) char(l) 

write (10,200) char (1) (1 : 4) , char ( 1) (5 : 6) , char(l)(7:8) 
, char (1) (9:10) ,char(l) (11 : 13) , char (1) (14:15) 

, char (1) (16:17) , char (1) (18:21) 
format (a21) 

format (lx, a4, 5x, 3 (3x, a2) , 5x, 2x, a3, 2 (3x, a2) , 5x, lx, a4) 
go to 10 

print *, 'Finished' 

stop 

end 
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10 


* 

■k 


100 

200 


•500 


character*29 char (1) 
read ( 9, 100 , end=500) char(l) 

write (10,200) char (1 ) ( 1 : 6) , char (1) (8 : 11) , char (1) (12 : 15) 

, char (1) (16:18), char<l) (19 : 21) , char (1 ) (22:23) 

, chard) (24:25) ,char(l) (2 6 : 27) , char (1 ) (28:29) 
format (a29) 

format (lx, a6, 5x, 2 (lx, a4) , 5x, 2 (2x, a3) , 2 (3x, a2) , 5x, 2 (3x, a2) ) 
go to 10 

print *, 'Finished' 

stop 

end 
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